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Abstract 

Calculation aspects of holonomic quantum computer (HQC) are considered. 
Wilczek-Zee potential denning the set of quantum calculations for HQC is ex- 
| plicitly evaluated. Principal possibility of realization of the logical gates for this 

case is discussed. 

m 
o 

(Nj | The conception of quantum computer (QC) and quantum calculation developed in 

80-th [I] [2] were found to be fruitful both for computer science and mathematics as 
well as for physics Although a device being able to perform quantum calculations 
is far away from practical realization, there is a number of theoretical proposals of such 
Qh! a construct [I] [IE]. Recently holonomic realization of QC was proposed [E|,[IB]- It is 

based on the notion of non-Abelian Berry phaseJHU-To perform quantum computations 
in this approach one needs a parametrically driven quantum system described by the 
Hamiltonian H(X) with adiabatically evolving parameters X(t) = (Ai(£), . . . , Ajv(£))- 
Adiabaticity means that Q <C ^min where Q is the characteristic frequency in the 
Fourier spectrum of X(t) and uj m i n is the minimal transition frequency in the spectrum 
of H(X). If the spectrum is degenerate then the cyclic evolution of N parameters A^(t) 
results in a unitary transformation of each eigenspace of H. Such transformations can 
be treated as computations on the eigenspace representing a part of the qubit space of 
HQC. These computations are shown to form a universal set of quantum gates for a 
specific relevant model j20]- For instance, one should choose of independent 2-loops C 
in the control manifold M = {all possible A} and then the non-Abelian Berry phases 
U (C) corresponding to each loop represent all basic gates. To produce such gates one 
should have a possibility to control i.e. for each closed adiabatic curve C in M one is to 
know the explicit form of U(C). To know that it is necessary to be aware of an explicit 
expression of non-Abelian Wilczek-Zee potential A which determines U(C) completely. 
A more or less universal method of explicit calculation of A is to our knowledge up to 
now absent. For an early overview on this subject see [21]. A variant of Berry phase 
evaluation for a finite level system relevant to HQC was proposed in ^H] where a specific 
parametrization of SU(n) is used. In this paper we present some less general but to 
our mind some more effective method of evaluation of non-Abelian Berry phase which 
can appear to be useful for HQC. For this purpose we restrict ourselves with systems 



*National Centre of High Energy Physics, Bogdanovich Str.153, Minsk 220040, Belarus 
t Institute of Physics National Academy of Sciences Skoryna av.68, Minsk 220072, Belarus 



1 



with finite number energy levels described by isospectral Hamiltonians, i.e. by ones 
whose spectrum does not depend on time via adiabatical evolution of their parameters. 
At first we demonstrate the method for the Abelian case, then derive the result for 
the non- Abelian one and in conclusion some possibilities of physical realization of the 
basic gates. 

In veew of above made restrictions the Hamiltonian we consider can be present in 
the form 

H{t) = V{t)H V\t), (1) 

where the n x n matrix Hq does not depend on time. As the qubit space must be 
large enough then the state space of H must be so as well. In this case derivation of A 
becomes difficult even if there are few independently evolving parameters that is the 
case for a typical experimental situation. 

We derive an explicit closed expression for A in terms of Hij(t) that have a direct 
physical sense. Let us consider the Hamiltonian H(X) with Aa(£) evolving adiabatically 
thus Born-Fock theorem [22] is valid: 

H{\)<p?(\) = Ei(\)<p?{\), z = l,...,s, a i = l,...^. (2) 

Here di is the degeneration rate and J2i=i d-i = n. It should be emphasized that 
the finite-matrix form of the Hamiltonian H(X) may be not only the result of an 
approximation, i.e. one neglects all energy levels whose number is greater than n but 
also a consequence of the symmetry of the system under consideration. This means if 
the Hamiltonian H we start with possesses a Liean symmetry, i.e. it can be represented 
in the form H = ca(A) Ta, where T 'a are generators of a Lie algebra, then the state 
space of the system is split in a set of components irreducible with respect to the 
algebra and for each finite-dimensional component we come to the finite-matrix form. 
Cyclic evolution of A results ^H] i n a monodromy transformation of ip = ((pj, . . . , y^) : 

cp(T) = U{C)<p{0), (3) 
where t parametrizes the curve C so that A(T) = A(0). 

For an illustration we start with Abelian Berry phase that was firstly predicted by 
Berry [22] • in this case all di — 1 and instead of (J3J) the relation 

<^ m (T)=e^( c V m (0), (4) 

where 

7m(C) = J A m , A m = i (v? m \d(p m ) . (5) 
c 

As the spectrum of the Hamiltonian is finite, the ket \ip m ) is a unit vector m in C n , 
so A m is 

i 

A m = - (m*dm — mdm*). (6) 
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As the evolution is adiabatic, the spectrum of H(t) remains always nondegenerate if it 
was so at the initial time. Then there is always a nonzero main minor of H — E m which 
we assume to consist always of the first n — 1 lines and columns of H — E m . Denoting 
the matrix consisting of the first n — 1 lines and columns of H by H± we come to the 
condition 

det(H ± -E m ) ^0 (7) 
Making use of this condition one can represent n in tyhe uniform coordinates 

m = — , 
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V 1 + 1& 

and express £ m in terms of for 1 < i, j < n — 1 and E m : 

£ m — {H± — E m ) 1 h, hi = —H in , (8) 
where h is a vector in C n_1 but not in C n . Thus we have for A m 

A m = ~ 1 i ic 12 ' 
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where £ m is completely determined by (JSJ). Note that the result obtained is purely 
geometrical because it can be expressed of the Kahlerian potential 

F = log(l + |z| 2 ) 

which determines all the geometrical properties of the state space 

SU(n)/U(l) x ... x U(l) . 

V v ' 

Ti—1 times 

Here z is a vector in consisting of n(n — l)/2 independent components of all £ m . 
One more observation is that for the case n = 2 corresponding to spin in the magnetic 
field © gives well known result [21 



2 J 1 + |e| 2 2 V ; ' 

c 

where ± labels spim-up and spin-down states, £ is thought of as a stereographic pro- 
jection of a point on the unit sphere S 2 and Q(C) is the solid angle corresponding to 
the given closed contour on the sphere. 

Now we proceed with a more general case of degenerate spectrum. Quantum compu- 
tation for this case generated by an adiabatic loop in the control manifold is determined 
by © where U(C) is presented by a P-ordered exponent 



U(C) =Vexp^ A, 



{A m ) ab = i(m b \dm a ) . (10) 



The set of eigenvectors £ ma , a = 1, d m must obey the equation 

{Hj_ — Em) £ma = ^ C a- 



Here the matrix H^f m ^ is constructed from the first n — d m lines and columns of H, c a 
are arbitrary d m - dimensional vectors and h is the following (n — d m ) x <i m -matrix: 



d m +l 



\H n —d m ,n—d m +l ■■■ H n _,i mn J 

Of course it has sense only if the condition 

det(H^\t)-E m )^0 



(12) 



is valid along the evolution process. The set of vectors £ ma must be orthogonalized 
by the standard Gram algorithm and after that we get the orthonormal set of the 
eigenvectors z a (here and below we has omitted the index m) in the form 
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/ 



det T a _i 



x 



A 



rvi ; 



(13) 



where x b = (£ fe , c b ) and c a is chosen to be the standard orthogonal set c a = (0... 1 ...0). 
The matrices r a are determined by 



i + ztz a , 



(14) 



where the (n — d m ) x a-matrix Z a consists of a first lines of the (n — d m ) x rf m -matrix 
Z = (H^f™^ — E m )~ x h. Using ()14|) and ()13j) we come to the final expression for the 
matrix-valued 1-form A: 



2 det(l + Z\_ x Z a _ x ) det(l + Z\_ x Z h _ 



1 < % < a, 1 < j < b, (15) 



where 



-pi -p*J 
1 a L b ' 



^ab 



min(a,b) 

(£ j \d1m(g$)\t i )+ rfIm (^), 



i=l 



and is the cofactor of ^ in T a . Note that the change of our basis c a by c' a = U a b(X) c b 
leads to a standard gauge transformation of A 



A' = UAU ] + i{dU)U^ 
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The formula (|15|) is the desired expression of A in terms of the matrix elements of the 
Hamiltonian. It is correct if condition (|12|) is valid. It is not nevertheless a principal 
restriction because d m does not depend on time due to adiabaticity of the evolution 
and there is always at least one nonzero n — d m -order minor of H . Then, if the minor 
we choose vanishes somewhere on the loop C one can always take local coordinates 
such that the techniques considered is applicable on each segment of C. 

One more thing to discuss here are quantum gates and their possible physical re- 
alization. As the Hamiltonian is an element of a Lie algebra which we assume to 
be semisimple, generic computations are of course elements of the corresponding Lie 
group. We represent them in the form 



where are real and ( a are complex parameters, k is the dimension of the Cartan 
subalgebra, i labels all it's linearly independent elements, the product is taking over all 
positive roots of the algebra, and the standard notations for Cartan-Weyl basis have 
been used. As the action of the first exponent in (|16|) results in a trivial phase factor, 
a generic computation is determined by the product over all roots. Thus we come to 
the conclusion that the elementary factors 



provide the basis of quantum gates for the model. Their number is e.g. for the first 
fundamental series A n n(n — l)/2 as it should be. An idea of physical realization of 
the gates appears if one takes into account that for A n E a can be realized by means 
of ordinary bosonic creation and annihilation operators, namely E a = a^cij for some 
1 < h 3 ' < n. Then E a represents nothing but two-mode squeezing operator. Thus the 
model considered can be applied to optical HQC with n laser beams (the case n=2 is 
considered in j2U|) and the logical gates U a are just two-qubit transformations realized 
by transformation of two laser beams. 

The method presented here enables one to build in principal any computation for 
HQC described by a Hamiltonian with a stationary spectrum in terms of experimentally 
measured values, namely the matrix elements of the Hamiltonian. The method depends 
weakly on the dimension of the qubit space which other models based on various 
parametrizations of evolution operators of the system are very sensitive to. Application 
of this method to a concrete physical model will be discussed elsewhere. 
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